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I. INTRODUCTION 


For the Army to fight effectively in a resource scarce 
environment, the quantitative decision making techniques of 
Operations research are important skills for Army staff 
officers. Staff officers are expected to be able +o put 
humbers in their estimates when briefing commanders. They 
are expected to be able +o measur2 and evaluate complex 
operations and subordinate units. They are expected to be 
frugal managers of time and money. And above all, staff 
officers must be able to apply sound, quantified reasoning 
in planning how to win the air-land battle. 

The use of hand-held programmable calculators by Army 
staff officers has the potential for improving the use of 
quantitative decision making technigues throughout the Army. 
Faster and more accurate than paper and pencil, the calcula- 
tor is less expensive and more portable than larger comput- 
ers. Even when compared to the latest micre-computer 
systems or to portable terminals used for diszributed data 
processing, the hand-held programmable calculator offers 
advantages in cost, reliability, power consumption and emis- 


Sion of electromagnetic radiation. Hand-held programmable 





calculators have already been successfully used by soldiers 
moe che field for applications in artillery Eire direction, 
surveying, and navigation. In addition, large numbers of 
Acmy officers own their own pocket calculators and routinely 
use them for staff planning and reporting functions. 

In January of 1987 the U. S. Army Command and General 
Staff College at Fort Leavenworth, Kansas selected a pro- 
grammable calculator for the Combined Arms and Services 
Staff School (CAS 3 .) USing both resident and non-resident 
instruction, this course 1s deSigned to teach ali Army cap- 
tains staff techniques and procedures. AS a Significant 
part of the curriculum, the students are introduced to sub- 
jects such as statistics and regression, decision theory, 
combat modeling and linear programming. Considering the 
large number of officers projected to attend this course in 
future years, this course represents the most widespread 
training in operations research techniques ever attempted by 
the Army. The decision to provid? a sophisticated calcula- 
tor to these students on an experimental basis was made for 
two fundamental reasons. First, the availability of a cal- 
culator with immediate field utility should motivate the 
Student to apply the quantitative techniques as compared <*o 


the student who would be forced to do all calculations by 





hand. Second, the power of the calculator permits classroon 
discussion cf techniques such as linear programming and 
regression which are very difficult and time consuming to 
perform manually. 

This thesis documents the author's work to support the 
use of a calculator in the Combined Arms and Services Staff 
School. Initially, the intent was to produce a series of 
lesson materials incorporating the use of the calculator on 
a series of operations research topics which have immediate 
application for the Army division level staff officer. 
Instead, the work accomplished focused on the design and 
construction of a system to make the programming and testing 
of calculator programs easier and more efficient. Except 
for the introduction, this thesis is written for the person 
wishing to implement the pregramming support systen 
described. The implementor must have a detailed knowledge 
of the instruction set and programming characteristics of 
the HPYiCV calculator as described in Wickes [Ref. 1: pp. 
6-20]. For the eventual user of the system, as compared =o 
the implementor, the system itself provides on-line dccumen- 
tation on how to use the system and what commands and 
options are available. Figure 1 Shows the command menu dis- 


played on the terminal screen by this interactive progran; 
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Figure 2 gives a more detailed explanation of each of the 
commands; and Figure 3 displays the on line introductory 
material that is provided to new users of the system. For 
the user, a knowledge of the information contained in the 
calculator owner's handbook [Ref. 2] is sufficient to begin 
writing calculator programs using the support system 
described. 

The calculator selected by the Command and General Staff 
College, the Hewlett-Packard HP4YICV, typifies the state of 
the art in off-the-shelf calculator technology. While not 
Without disadvantages, this calculator was selected because 
of its power and features which make it easier for Arny 
staff officers to use. First and most important of these 
features is the ability of the calculator tc manipulate 
alphabetic characters in addition to numeric data. The cal- 
culator can display the name of a variable when input data 
is required or label output when the calculation is con- 
pleted. With this feature, the calculator helps the user 
know what data to input or what action to take next. It 
also helps alleviate the need for constant reference to 
printed instructions which are difficult to use under field 
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A second important feature-is the multiplicity of means 
by which programs can be entered into the calculator. Mag- 
hetic cards, read only memory, and optical bar code are ail 
available and each has advantages depending on the situa- 
tion. Fer the long term, read only memory offers the abil- 
ity to retain very large programs (in excess of 8000 bytes) 
and the Simplest and most reliable means of entering pro- 
grams into the calculator under field conditions. For the 
Short term, optical bar code offers the least expensive 
method of reproducing and distributing calculator software 
that has not been subject to extensive field testing. In 
addition, as shown in this thesis, the optical bar code can 
provide an important link between a main-frame computer and 
the hand-held calculator. 

A third important feature of the HP41CV is its rela- 
tively large memory capacity as compared to programmable 
calculators such as the Texas Instruments TI-59. A large 
amount of memory permits the solution of larger, often more 
realistic problems than could previously be solved ona 
hand-held device. A demonstration program given in this 
thesis for linear programning is an example of an applica- 
tion where the full memory capability of the HPYICV is 


required =o be able tc solve realistic problams. 
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To take advantage of the calculator's unequalled economy 
and portability, the operations research analyst is chal- 
lenged to overcome its limits of speed and memory capacity. 
Rmempreparation of calculator software is as difficuit, if 
not more so, than the preparation of software for larger 
computers. To accomplish the most possible with the hand- 
held device, the calculator programmer is often forced to 
write programs which are very difficult to comprehend when 
examined by other programmers. As Dahl, Dijkstra and Hoare 
exec. 3: pp. 1-10} point out there are limits to human con- 
petence which interfere with «he programming process. [In 
the past, with less mature calculators which constrained the 
typical program to a few hundred program steps, these limits 
to human competence were neither as apparent nor as economi- 
cally important as they are with the HP41CV. Accordingly, 
it is not envisioned that the average Army officer who uses 
the HPYICY on real world problems which push the calculator 
momecne limits of its capability would write their own pro- 
grams. In particular, it was never intended that the stu- 
dents in the Combined Arms and Services Staff School would 
be taught calculator programming. It is a tribute to the 
power of the device and the quality of the calculator soft- 


ware when a relatively inexperienced user can run complex 
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programs using little more than the digit entry kevs and che 
crun-stop key on the calculator. This does not mean that the 
user must not have a clear understanding of his problem or 
the solution technique, but rather it means that the calcu- 
lator should not require programming skill or extensive 
meaiieng prior to application. 

The growing complexity of caiculator programs described 
above and the realization that calculator programs for Army 
field usé are not programmed in the field, suggest the need 
for a system to Support the development, distribution and 
maintenance of calculator software. An operations research 
analyst or other professional programmer must be able to 
more efficiently prepare calculator programs than by keying 
them into the hand-held device. By preparing the programs 
initially on a larger computer, such as the IBM 3033, the 
programmer can use the speed and storage capability of the 
larger machine to great advantage. In addition, the avail- 
ability ef a full-screen video text editor speeds the pro- 
cess of program revision and maintenance. By providing a 
capability to integrate comments directly into the source 
code on the larger computer, program documentation is more 
easily provided. Essentially the idea is that a programmer 


WOuld write the calculator program using a terminal 
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connected to a large computer. After the calculator program 
is entered into the large computer, a compiler program 
running on the large computer would check the calculator 
program for errors and convert the mnemonic instructions 
into the "key codes" which ar2 the numeric instructions 
actually executed by the calculator. Then an emulator 
program running on the large computer would take «he numeric 
instructions from the compiler and execute the program--in 
effect making the large computer produce the same effects as 
the calculator only much faster and more efficiently for the 
pregqrammer. Finally, when the program has been written and 
tested on the large computer, optical bar code is produced 
Which allows for the economical distribution and use of the 
program in the field. T9 encourage the calculator progran- 
mer to use the system described, this process should occur 
in an interactive programming environment in which the user 
can move from one step to another by issuing simple commands 
such as those listed and described in Figure 2 and receive 
help or cn line instruction whenever desired. Under this 
proposed system, the advantages of both the larger computer 
and the hand-held calculator are used appropriately ina 


mutually supporting manner. This thesis presents two of the 


components of this proposed system. First, an IBM EXEC II 
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program is given which provides an interactive programming 
environment for users operating under IBM's Conversational 
Monitor System (CMS.) A short discussion of the design of 
this program and a complete copy of the source code is 
contained in Appendix Dto this thesis. Secondly, a cross 
compiler written in IBM standard FORTRAN IV is provided for 
translating calculator mnemonic inscructions into the key 
codes necessary for use by the emulator and also for the 
Production of optical bar code. The term cross compiler 
refers to the fact that the program runs on one machine (the 
larger computer) but compiles programs for another machine 
(the calculator.) A discussion of the design of this 
program and a complete copy of the source code is contained 
in Appendix D. To make the program easier to understand and 
adapt to new requirements, it is modularized into 24 
subroutines and is heavily commented. 

To illustrate the use of the system, two of the six 
example programs originally planned are provided in this 
shesis. Revised plans now call for the remaining four exan- 
ple programs to be issued at a later date as Naval Postgrad- 
uate School technical reports. Because the reasons for the 
delay constitute some of the most important lessons learned 


from this thesis research, Chapter 2 documents the process 
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weecioeae ~e€Chnical discussion of the factors involved. The 
major conclusions described in Chapter 2 are the need for a 
@emorttized list of Criteria with which to evaluate calcula- 
tor programs and the need for mors structure in the progran- 
ming process. Chapter 2 is technically oriented and assumes 
+he reader is familiar with the concepts of structured 
programming. 

Each of the calculator program examples is described in 


a separate appendix in which the documentation listed in 


1. Program Description 


2. Sample Problem 


4. Source Code Listing with Comments 


| 

| 

| 

| 

| 

| 3. User's Guide 
| 

| 

| 5. Bar Code 
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Figure 4: Components of Program Documentation 


Figure 4 is provided. The first example on single variable 
Statistics is documented in Appendix A and uses the calcula- 
tor in an area where calculators have long been used, but 


does so in a way that shows the unique capabilities of the 
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HPYI1CV. <A second example on linear programming is docu- 
mented in Appendix B and illustrates an area where calcula- 
tors have not received widespread application. Most 
calculator linear programs which have been published to date 
have been either incomplete algorithms or have been limited 
to very Simple problems. 

A third example, which by its nature does not conform to 
the documentation standards outlined above, describes a set 
Semicality toutines which could be distributed in read only 
memory. Programs for read only memory have different char- 
acteristics from other calculator programs and Appendix C is 


provided to illustrate some of these differences. 
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A. CHAPTER OVERVIEW 

This Chapter examines calculator programming within the 
context of the author's experience in preparing HP4YICV pro- 
grams in support of the Combined Arms and Services Staff 
School. With the advanced capabilities and features of the 
HP4Y1CV, it was hoped that a complete package of software 
could be prepared quickly. To document why this did not 
occur, this chapter will examine strengths and weaknesses of 
the calculator in relationship to a collection of techniques 
referred to in computer science as structured programming. 
For the reader unfamiliar with this term, the previously 
emredq WOLK by Dahl, Dijkstra, and Hoare [Ref. 3]j is 
recommended. This chapter is technically oriented and does 
assume familiaritiy with structured programming concepts. 

When programming calculator programs for personal use, 


most programmers, including the author, do not find the task 


| 


ct 


Geectcult. Program@ang 2 hand-held calculator with the 
capabilities and features of the HPY1ICV can be a rewarding 


experience. It is rewarding +o master the algorithm of an 


Operations research technique on a hand-held device. The 


23 





educational value in programming the calculator has been 
recognized by many educators, including Hamming [Ref. 4: 
Gemee-5 | and Weir [8ef. 5: pp. xli=-xiiij. Providing a pro- 
gram for general distribution which makes optinum use of the 
calculator is quite a different situation. It was the 
author's experience that programs, which gave correct 
answers when used by che author, often had to be completely 
re-written several times before being acceptable. This 
problem became more acut? as the size of the programs grew 
beyond 400 program steps, for at that size it became 
increasingly difficult *to modify programs without affecting 
the total design. The major conclusions described in this 
chapter are the need for a prioritized list of criteria with 
which to evaluate calculator programs and the need for more 


Structure in the programming process. 


B. STRUCTURED PROGRAMMING WITH THE HPUICY 
1. The Need for structure 
To increase the efficiency of the programming pro- 
cess, a collection of techniques known as structured pro- 
gramming has received widespread attention in the computer 
science community. While there is no one definition of 


structured programming, it does require three essential 


characteristics. First, there must be a logical structure 
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+o the program which reflects the nature of the problem to 
be solved and any constraints imposed upon th2 soluticn. 
Second, the systematic process of stepwise refinement is 
used to limit the complexity of program segments. Third, 
the programming language must reflect the logical structure 
of the program and assist in stepwise refinement. These 
three characteristics represent not so much a detailed réc- 
ipe for program development as they do a philosophy of how 
programs can be more efficiently written. It was with this 
philosophy in mind, that a calculator programming support 
system was proposed which could take into account the 
strengths and weaknesses of the calculator; balance the 
Structured programming philosophy with the other criteria 
listed below; and thereby solve the problems encountered in 
Writing calculator software for the Compined Arms and 
Services Staff School. 
2. Fundamental Limitations of Calculators 

Writing programs to solve complex problems on a 
hand-held calculator is difficult both because of inherent 
limitations in the calculation speed and memory capacity of 
the machine and also the inability of the calculator's 
native programming language to directly support structured 


programming constructs. In many respects, the task is 
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similar to writing assembly level language programs fer 
larger computers. Calculator programming features a 
powerful instruction set including advanced mathematical 
functions but lacks any ability to refer to variables by 
name instead of storage address. Like assembly language, 
the calculator's programming language consists of short 
mnemonic instructions typically followed by the storage 
iecation of the data to which the operation is to be 
applied. While a large amount of computer programming is 
Still done in assembly language, it is generally accepted 
that programming in a higher level language such as FORTRAN 
is preferable. Programs written in an assembly language 
take more time to write and are not as easily changed as 
higher level language programs. Also, because they depend 
on the instruction set of a particular machine, they can not 
be easily transfered from one computer to another. These 
same disadvantages apply to calculator programming. In 
addition, because the hand-held device does not have the 
speed and memory capability of the larger machine, the 
calculator programmer must be even more mindful of the need 
to optimize his program to save program steps and execution 


time. 
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3. Modular Design 

The HP4YICV supports structured program@ing as well 
or better than any other hand-held calculator. As described 
in the owner's manual [{Ref. 2: pp. 177-196], the machine 
primitive instruction XEQ encourages the construction of 
modular programs using calculator subroutines. Each subrou- 
tine can be a self-contained unit capable of being written 
and tested independently and used by multiple programs. 
This modularity is most strongly encouraged when routines in 
read only memory are used, for then the application progranm- 
mer can significantly reduce the number of program steps in 
his own program. This modularity, however, is rot complete, 
Since all variables are globally referenced and can be 
changed daliberately or inadvertently by any subroutine., 
This problem 1S no more apparent than with the use of read 
only memory, Since one of the most limiting factors in using 
the read only memory programs as subroutines is conflict in 
the use cof common registers. Also, unlike the modularity 
required in truly structured programs, there is no 
restriction limiting a subroutine to a single entry anda 
Single exit point. In structured programs, such limits on 
entry and exit serve to define the fundamental building 


blocks by which stepwise refinement is made possible. With 
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the calculator, however, multiple entry and exit points are 
most useful for allowing a common routine to handle a 
@mplicity of problem conditions. In this thesis, for 
example, programs are given for which two standard entry 
points are provided. One entry point uses an alpha-numeric 
label and an audio prompt to speed data entry, while a 
second entry point uses the alpha-numeric label but 
Suppresses the audio tone. After data entry, the value 
entered is displayed, and the user is required to verify the 
accuracy of the data entered. By using the same subroutine 
with different entry points, memory space is saved overall 
at the sacrifice of the structured programming philosophy. 
4. Control of Program Plow 

A basic deficiency prohibiting the HP4YI1CV from 
directly supporting structured programming is the way in 
Which program flow is controlled. Programming languages 
which support structured programming typically have instruc- 
tion constructs such as WHILE--~ENDWHILE, REPEAT--UNTIL, or 
LOOP=-QUIT~-~ENDLOOP which make programming loops clear and 
concise. Constructs such as IP--THEN--ELSEIF<--ELSE--ENDIF 
and the CASE statement make the evaluation of conditional 
expressions efficient and relatively error free. Also, 


Structured programming languages typically discourage the 
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use of GOTO unconditional transfers because they lead to 
confusing code. In contrast, the HP4Y1ICV programmer must 
write his own looping constructs and his own conditional 
evaluation constructs using machine primitive instructions 
which somewhat obscure the program's basic objective and 
Mmmeeumot CONtLOL, In addition, it is dafficult to avoid dis- 
turbing pending operations in the stack registers when a 
conditional statement must be evaluated. As can be seen by 
the short program shown in Figure 1, the notation of the 
programming language does not permit structured program 
iLOW . 
5. Clarification of Program Structure 

Because no calculator, including the HP4Y1ICV, sup- 
ports named variables, the use of comments as an integral 
part of the calculator program is vital if the logical 
structure of the program is to be made clear as required by 
structured programming. Comments should provide the vari- 
able names when storing and recalling data; they should pro- 
vide clarification of program flow; and they should mark 
subroutine boundaries and entry and @xit points to make it 
easier to identify segments of the program. With the 
HP4YICV's stack oriented architecture, it is also frequently 


useful to display the names of the contents of each of «he 
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fragment will sum the data values stored in memory 


locations 1 through n. 


| 
| Given che number n in the x-register, this program 
| 


instruction Commen® 
LBL "Svs To execute press "XEQ SUM". 
ves 
y. 
1 
+ Establishes a loop counter. 
0 Clears x and pushes loop 
LBL 00 Gountos pa neo ry < 
RCL IND Y Recall the next data value. 
+ Accumulate the sun. 
RSG nL Increment the loop counter. 
GTO 00 If more data remains, branch; 
RTN else, quit and display sun. 


Figure 5: 
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Example Program to Add n Numbers 
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stack registers. In Appendix C on common subroutines with 
read only memory application, a shell sort [Ref. 6: pp. 
84-95] routine is given which employs the technique of using 
comments to display the names of the variables on the stack 
register. 


Ge ba & 


Types and Indirect Addressing 

Calculator programs represent more than a sequence 
of keystrokes; they also represent the manipulation and 
transformation of data. For maximum efficiency, the manipu- 
lation of data should be structured so as to prevent common 
programming errors. For this reason, most computer lan- 
guages which directly support structured programming enforce 
data type correspondence between data and operations. Fre- 
quently the formal declaration and initialization of vari- 
ables is also reguired. The HPYICV handles two types of 
data--real numbers and alphanumeric characters. While no 
formal declaration of variables is required, «type checking 
is done automatically and is transparent to the user. Any 
attempt to perform an arithmetic operation on alpha-numeric 
data will result in the message "ALPHA DATA" and the program 
wee) halt. 

Because there is no formal declaration of variables, 


the programmer writing programs for the HP41ICV must use 
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axtreme caution in managing his data set and insuring that 
the numbers stored and recalled by the calculatcr progran 
are in fact the data elements desired. A typical example of 
an improper data reference occurs when a pregram is using 
indirect addressing and attempts to store or recall data 
from a non-existent data register. This programming error 
is so common that a special error message "NONEXISTENT" is 
provided by the calculator when this error is detected. 
Indirect addressing is an important feature which gives the 
calculator a considerable amount of power and flexibility, 
but also represents an additional responsibility for the 
programmer to explicitly control. On the HP41CV all indi- 
rect addressing calculations must be specifically provided 
by the application program--there are no vector or array 
data types such as usually found with higher level lan- 
guages. In an attempt to make indirect addressing more 
transparent to the programmer, an experimental subroutine 
WaS prepared to recall an arbitrary element of a matrix 
stored as a two dimensional array. This subroutine, which 
is shown in Figure 2, was used in a simultaneous differen- 
tial equation combat model and the results evaluated. It 
accomplished the task, but slowed the execution of the pro- 


gram considerably (resulting in an overhead of 10.5 seconds 
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of extra execution time for every 100 Subroutine calls) and 
did not significantly improve the size or legibility of the 
application program. Accordingly, this technique is not 

recommended and indirect addressing remains a task that must 


be treated explicitly by the application programmer. 


C. ADDITIONAL CRITERIA FOR PROGRAM EVALUATION 

Calculator programming in many respects resembles a 
multi-criteria decision problem. On the surface the 
criteria for program effectiveness are quite straight 
forward--the program must yield the correct answer, run 
quickly, require the fewest possible memory registers and be 
user friendly. Unfortunately, these objectives often 
conflict and can not always be simultaneously achieved. In 
particular, the principles of structured programming are 
often in conflict with the desire to reduce the size of 
programs and increase their execution speed. It is also 
crue that the objectives of structured programming concern 
the process of writing programs, whereas the additional 
@eat6ria listed concern the final program product itself and 
are therefore logically considered separately. Attempting 
to achieve all criteria at once can lead to failure, and 
some tradeoffs must be considered to evaluate programs and 


guide program development. The following criteria represent 
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Figure 6: 


Entry to this routine assumes the x register contains 
the column number and the y register contains «the row 
number. The base address must be stored in ROY and 


the dimension of the matrix must be stored in ROS. 


LBL "RCLM 
RCL O4 

+ 

X<>Y 

1 

ReLse.> 

~ 

+ 

ReLe END f 
RE 

END 








(BASE ADDRESS REGISTER 
(ADD BAS& TO COLUMN NUMBER 


(RECALL THE ROW NUMBER 


(DEMRENSZLONDOFP DTHE MATRIX 


(ADDRESS IS NOW IN X REG 


(RECALL THE DATA DESIRED 
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Program to Recall an Element of a Matrix 
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"lessons learned" in develofing application prodrams as 
examples for this thesis. 
1. User Friendliness 

User friendly programs consider the application 
environment and do not task th2 user to be all Knowing or 
Matnout Error in entering data. While individuals differ 
greatly with experience, the average user will make frequent 
errors in entering data with the hand-held calculator's 
small keyboard. In talking with officers who had used the 
Ge=59 Calculator in the field for fire direction, it was 
discovered that most preferred to use the printer with the 
calculator because it allowed data to be checked after 
entry. This was in spite of the fact that the printer and 
Calculator combination is more costly, less portable and 
less suitable for use in the field than the calculator 
alone. In short, user friendliness was more important than 
these other criteria. For this reason, it should be manda- 
tory that any calculator programs intended for Army use in 
the field must allow the verification of data after entry. 
Because the use of the printer obviates many of the advan- 
tages of the hand-held calculator, the printer should not be 
reguired for this verification. Ine of the considerable 


advantages of the HPYICV is that the large amount of prograan 
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memory makes it possible +o store the input values and 
perform this verification. However, programs written with 
this criteria in mind may not appear be most efficient +o 
the casual observer. 

Another important aspect of user friendliness is 
limiting the complexity of the calculator and the actions 
required =o get results. The typical Army officer has lit- 
tle appreciation for the multitude of scientific and mathe- 
matical functions labeling the keys of the HPYICV. Yet the 
common programming practice of using the toc two rows of 
calculator keys to indicate the identity of a variable 
elther upon input or output increases confusion over the use 
o£ the function keys. This works as follows: When local 
alphabetic labels are used in a program to represent entry 
points by which a user indicates the identity of an input 
Variable or requests a particular output variable, then «he 
first two rows of keys on the HPYICV become subroutine exe- 
cution keys pointing to these local labels when the calcuia- 
tor is in user mode. This feature was very important on the 
HP67 and TI-59 where the lack of alpha-numeric capability 
reguired this method of program execution in order to most 
easily determine the identity of the input or output value, 


but it is less important on the HP4ICV. It is almost always 
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true that a program which requires the use of local labels 
is harder to use, and requires more [frequent reference to 
the user instructions than a program which uses only the 
run-stop key and properly prompts the user and labels output 
values. 
2. Execution Speed 

The second most important criteria for a calculator 
program is that it must yield results relatively quickly. 
In preparing example programs for this thesis, this point 
became very clear when testing two particular programs. One 
program, a Simultaneous differential equation combat model, 
required in excess of 150 data values in order to yield 
results. It should be noted that it was only with the 
introduction of the HP4&ICV that it became feasible to con- 
Sider such large problems on a hand-held device. To accomno- 
date the size of the model, the program was written so as to 
economize on program steps at the expense of increased exe- 
cution time. It became immediately obvious upon initial 
testing that this had been the wrong priority--for users of 
the program were not impressed with either the use of the 
calculator or the utility of the combat model. If such user 
acceptance is not present, then the calculator program will 


remain unused, no matter how elegant the design to conserve 
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memory. In contrast, the linear programming example given 
in Appendix B was written so as to emphasize speed even if 
it meant including code redundancy. This program has been 
well received in part Beis it is so much faster than 
paper and pencil methods. 

The easiest and most effective technique that is 
useful in increasing speed is to decrease the number cf pro- 
gram steps that the calculator must process inside program 
locps. For example, if two different program options 
require similar but slightly different actions within a pro- 
gram loop, it is tempting to insert a program flag check and 
branching instructions within a loop so as to use the same 
loop for both conditions. But this means that the calcula- 
tor must ¢est the flag and branch inside the loop even 
though the program is probably shorter overall.! 

Instead, if the application permits, the memory capacity 
or the HP41CV can be used to best advantage by testing the 
flag once and then providing separate program loops for the 
two conditions. Again, this does not appear 2legant to the 
Casual observer, but it may result in a more successful 


program overall. This principle was discovered while 


1 Branching is required when the flag tests either set or 
Clear if more than one instruction iS required to account 
for the differences in the two conditions. 
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programming the single variable statistics program given in 
Appendix A. Initially, this program used a common loop for 
all data input and output operations, including reviewing 
the input data and making individual corrections. By 
providing a separate, somewhat redundant loop for data 
correction, the time required to input data points was 


reduced. 


D. A PROGRAMMING SUPPORT SYSTEM 

Considering the structured programming philosophy dis- 
cussed above in paragraph B and the additional criteria for 
evaluating programs listed in paragraph C, it becomes imme- 
diately obvious that programming with the calculator alone 
Will never meet even a majority of these objectives. [I+ 
must be recognized that the problem under consideration is 
not how the average person who owns a calculator should pro- 
ceed to program it for his own personal use, but rather how 
the Army can best provide the most cost-effective computa- 
tional resource for field use. For these reasons, a 
comprehensive programming support system is required. The 
programming support system outlined here will consider only 
the requirement for cost-effective preparation and 
Maintenance of the calculator programs and not the brceader 


issues of distribution and logistic support for the entire 
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calculator system to include hardware, training materials 
and printed references. 

1. A Cross Compiler and Bar Code Generator 

The first requirement for an operational support 

package is to free the programmer from the limitations of 
the hand-held calculator itself. Even with the printer and 
other peripherals, the calculator is no match for the larger 
Machine when large programs must be examined or edited. In 
addition, the calculator is not currently capable of produc- 
ing its own optical bar code as required for aconomic repro- 
duction and distribution of the software. Accordingly, a 
cross compiler for the HP41CV was listed as the first 
requirement of the programming support system. Such a cross 
compiler has been written and is the major outcome of this 
thesis effort. This cross compiler accepts an HP4YICV pro- 
gram writcen in the language of the calculator and returns 
the finished bar code as output. Any valid HP4Y1CV pregranm 
Will be processed without need for modification by the cross 
compiler. In addition to the basic language of the calcula- 
tor, the user is allowed to inject comments directly into 
the source code with the use of the left parenthesis as a 
comment indicator mark. The ability to make comments 


directly in the source code makes the calculator programs 


i 
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more legible and more easily modified at a later date or by 
another programmer. Often, well placed comments can make up 
for a lack of structure in the program itself as far as 
legibility and maintainability are concerned. Having the 
comments directly in the source code facilitates their use 
and helps insure that they are as up to date as the progran. 
For the average programmer, use of unmodified HP4YICY source 
code augmented with a comment indicator will represen+ the 
most common use of the cross compiler. The cross compiler 
is described in more detail in Appendix D including a 
complete listing of the source code. 

2. A Calculator Emulator : 

After the calculator source code has been processed 
by the cross compiler, a need exists to be able to run the 
program without the wait for the generation of bar code. [In 
addition, for the future development of read only memories 
for the calculator, an emulator program is required because 
the calculator itself can store only up to 2000 instructions 
in active random access memory. The read only memory can 
store up to four times this amount. Thus, the calculator by 
itself may not be capable of testing extremely large pro- 
grams or programs with large amounts of constant data also 


Stored in the read only memory. Although an 2mulator was 
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Moeewritten for this thesis, che design of the cross 
compiler reflects the need for such a program. For example, 
the cross compiler generates an intermediate array of deci- 
Mal integers which repesent the machine language of the 
HP41CV prior to conversion to binary. It was intended that 
these decimal integers could be used without modification or 
further translation within a FORTRAN computed goto state- 
Went. Thus, with the difficult translation, instruction 
parsing and syntax recognition already performed by the 
cross compiler routines, the emulator could consist of one 
large FORTRAN loop wherein a decimal integer was addressed 
in the instruction array by a program pointer variable. The 
integer is then immediately sent to a computed goto state- 
ment which would branch to the appropriate line of FORTRAN 
code which would simulate the referenced instruction, 
including updating the stack and the program pointer as 
appropriate. 
3. A Higher Level Language Compiler 

The final component in the calculator programming 
support system would be a program that would translate a 
higher level language such as PASCAL into HP41ICV language 
which could then be sent to the cross compiler for verifica- 


tion and generation of the bar code and intermediate 


42 





Savculator language listings. It is the higher level 
language compiler that would most directly make up for the 
weakness of the calculator in supporting structured progran- 
Ming. It would be able to increase the modularity of pro- 
grams, provide for named variables, make indirect addressing 
transparent and provide structured statements such as 
MOLLE--ENDWHILE and IF--THEN--ELSE. Again, the design of 
the cross compiler anticipates this requirement and provides 
a conSideraktle number of subroutines that would also be 
required by a higher level language compiler. These subrou- 
tines include a complete set of string functions for manipu- 
lating character data in FORTRAN and an instruction parser. 
Because it was envisioned that the higher level language 
compiler would also be able to process statements entered 
directly as HPYICV instructions, the cross compiler is con- 
structed so that the routine which compiles individual lines 
of HP41CV source code could be called as a subroutine by the 
higher level language compiler. Thus, all three major con- 
ponents of the proposed calculator programming support 


system would work together efficiently. 


43 





Pi KODUCTION: 

Calculating single variable statistics is one of the 
most frequently used applications of programmable calcula- 
tors. Army division level staff officers use single vari-~ 
able statistics to summarize and describe data for command 
briefings and periodic reports. The text by Mendenhall, 
Scheaffer and Wackerly [Ref. 7: pp. 3-13] 1s recommended as 
an introduction to the statistical measures calculated by 
the program given in this appendix. This program 


automatically calculates: 


e Mean and Median 


e Sample Standard Deviation 


e Sum of the Squared Deviations about the Mean 


e Coefficients of Skewness and Kurtosis 


e Minimum, Maximum and Range 


e Histogram Cell Frequencies 
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A single variable statistics program has been given as 
an example because of its immediate utility to the staff 
officer and to illustrate several features of the HP4YICV 
which make it a superior device for Army field use. The 
most important of the these features is alphanumeric rfrompt- 
ing for input data values. The program given in this 
appendix provides an alphanumeric prompt for every input and 
output value and requires only the digit entry kevs and 
run/stop key for data entry. Another important feature of 
the HP41cV used by this program is its large nemory 
Papacity. This program retains up to 219 data points in the 
calculator's memory to allow «he user to review the input 
data and make corrections during data entry. The large 
amount of memory allows the calcuiator to sort the data and 
calculate the order statistics including the mininua, 
maximum and median. Calculation of the median is a feature 
of this program which distinguishes it from other calculator 
statistics programs. In addition, without having to 
re-enter the data, the histogram may be calculated with a 


varying number of cells or a varying cell width. 


PROGRAM DESCRIPTION: 
The single variable statistics program has entry points 


for two different technigues of data input. The fastest 
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method, which provides both an alphanumeric prompt and an 
audio tone to speed data entry, may be called by execution 
of the program from entry point "STAT1." A slower method, 
which provides greater accuracy and suppresses the audio 
tone for classroom use, may be called by execution of the 
program from entry point "S11." When called from "S1," the 
program requires the verification of each data point after 
entry. The sequence of actions is as follows: 

1. The calculator displays an alphanumeric prompt. As 


an example, 'X1?% is the prompt for the first point. 


2. The user enters the data value with the digit entry 
keys and presses the rtun/stop key. 


3. The calculator displays the data entered with a label 
derived from the alphanumeric prompt. For example, 
alee ceo pacal caleulator response. § Thus 
display 1S prompting the user to verify the correct- 
nesS of the data displayed. 


4, If the value is correct, then the user simply presses 
the run/stop key and the calculator advances to the 
n2xt point. 

5S. If the value is erroneous, the user enters the correct 
value with the digit entry keys and then presses the 
run/stop key. Then the calculator will again repeat 
step 3 and ask the user to verify the data value. 

This process will continue until the user makes no 
modification to the data value. 


To run the program from either entry point the user may use 
the XEQ key, or asSign the entry point label ("STAT1" or 
"S1") to a key and execute it by pressing that key in the 
USER mode. Further instructions on running programs and 
Making key assignments are contained in the calculator 


owner's manual (Ref. 2: pp 114-116]. 
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Mca fone tO ties two initial entry points described 
above, several other alphabetic labels provids the user with 
functions that are called outside the normal sequence of 
program execution. Label "SR" provides the user with the 
capability to review the data stored in calculator memory, 
@ither before or after the data has been sorted. When used 
before the sort, the "SR" function is most useful in verify- 
ing the entire data set at one time. If used for this pur- 
pose, it should be called after all of the data has been 
entered and the mean of the data set is displayed with the 
"XYBAR" label. If flag 21, the printer enable flag, is set 
"'on" during this data review, then the calculator will stop 
as each point is displayed and the user may make corrections 
in the same manner as described above for the point-by-point 
verification associated with the "S1" entry point. When 
Mmeeag after the sort, the "SR" function is most useful for 
displaying the order statistics for the data set. If used 
for this purpose, it should be called after the histogram is 
output--when the "CMD" prompt is displayed. If the user 
presses run/stop after the "CMD" prompt, the order 
statistics will automatically be displayed. 

The design of the program, especially the data entry 


loop, reflects the need for calculation speed. Code 
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Bedundancy exists at several points in order £0 reduce «he 
need for extra flags, labels and goto statements which would 
slow execution during data entry. In spite of this need for 
speed, the summary totals needed for calculation of mean, 
standard deviation, skewness and kurtosis are accumulated 
during data entry. This is done so that these summary 
statistics are available with little or no wait following 
data entry. 

A complete listing of the program registers and flags 
used by this program is shown at the end of the program 


esting. 


SAMPLE PROBLEM: 

In order to establish a training standard for an obsta- 
cle courses, a division assistant $3 randomly selects 10 sol- 
diers and records the time it takes each to complete the 


course. The following times in minutes were recorded: 


Determine the summary statistics and cell frequencies 


hecessary to plot a histogram of this data. 
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SOLUTION: 


1. First, set the size of the calculator's data nemory 
large enough to retain the data values. This requires at 
least 16 registers plus 1 for each data point, or a total of 
26 in this example. Alternatively, the size of data memory 
may be set arbitrarily large, up to a maximum of 235 pro- 
vided the user has no other programs in the calculator he 


Wishes to retain. For this example press: 
XEQ ALPHA Sos ALPHA Zo 


2. To call the program, determine the appropriate method of 
data entry and select the corresponding entry point. 
Press: 


XFQ ALPHA STAT1 #£ALPHA (quick entry) 


XEQ ALPHA S 1 ALPHA (classroom use) 


gerne Calculator will respond with the prompt "N?" asking 


for the number of data points. Press: 
10 R/S 


feeeernie Calculator will vespond with the prompt "X11?" asking 


aa 


Bee the first data point. Press: 


2.1 R/S 
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Seen you called the program via "S1" the calculator will 
respond with "X1=2.100" asking for verification that the 
Mees t DpOUNtT 2S COLrrect. Lf not correct enter the correct 


value, else press run/stop. 


6. The calculator will continue in the same way as steps 4 
and 5 for the remaining data points until all the data has 
been entered. If at any time you discover that you have 


made an error in data entry for any point, pr2ss: 


X EQ ALPHA SC ALPHA 


The calculator will respond with the prompt "POINT?" asking 
for the number of the point in error. For example, if point 


number 5 were in error, you would then press: 


5 R/S 


Assuming you had just input a5 as the point in error, the 
calculator would then respond with the prompt "X5?" asking 
for the correct value of point 5. Respond with the correct 
value and press run/stop. The calculator will then go back 
to the place in the data entry sequence where it left off or 
1t will go to the calculation of the summary statistics if 


data entry was previously completed. 


50 





fameewnen data entry has been completed, the calculator will 


respond with the mean of the data sample labeled as follows: 


XBAR=2. 470 


At this point, you have the option of reviewing the entire 
data set or continuing to calculate the remainder of the 


statistics. To review the entire data set, press: 


XEQ ALPHA a ALPHA 


Hore that 1£ flag 21 1S set on (press SF 21), the calcula- 
tor will stop after each data point is displayed, permitting 
you to change any value simply by entering the new value and 


pressing run/stop. 


8. After the mean is displayed with the "XBAR" label, if 
yOu Simply press the run/stop key, the calculator will cal- 
culate the following statistics with the label shown: After 


each press R/S. 


Display Meaning 

SSQD=0.521 Sum of Squared Deviations 
About the Mean . 

SX=0.241 Sample Standard Deviation 

SKEW=0.170 Skewness 

KURTO= 2. 302 Kurtosis 


9. At this point the calculator will automatically sort the 


data. This may take from several seconds to several minutes 
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rH 


depending on the number of points in the data set. Afte 
the data set has been sorted, the calculator will display 
the median as follows: 


ae QO TO (Press R/S) 


2.800 
Two data values are displayed because when the number of 
data points is even, the median is not unique, but rather 
spans an interval from the one point listed above to the 
other. Many users may wish to simply take the middle of 
this interval as the median, but any point is technically 
correct in the interval. When the number of data points is 


odd, the median is unigue and only one value will be 


displayed by the calculator. 


10. After the median is displayed as described in step 9, 
the calculator will display the following statistics labeled 


as shown: 


Display Meaning 
MIN=2.100 Minimum Value 
MAX=2.900 Maximum value 
RNG=0. 800 Range 


mies At this point the calculator will respond with "CELL?" 


asking for the number of cells the user desires in the 
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histogram. If the number of cells is not Significant at 
mas point, the calculator will pick an appropriate number 
if the user simply presses run/stop. For this example, 


press: 


Kyo 


12. Next the calculator responds with "WIDTH" asking for 
the width of the cells. Simply press run/stop if you do not 
Wish to establish the width manually. Again, you may see 
the width the calculator will use by pressing the clear 
arrow key (Unless the width 1s an integer, you will also 
need to press FIX 3 to display the decimal properly if you 


wish to examine the width.) For this example, press: 


R72 


13. The calculator will now display the cell frequency 
counts as an integer count followed by the next ceil bound- 
ary. The leftmost cell boundary is set equal to the mininun 
value and is not explicitly output. If a data point falls 


exactly on a cell boundary, it is counted in the lett cell. 
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For this 2xample, the display will show: 


Display Meaning 

CNT=2 Iwo, ODSCEVat2oOnS i... 

between 2.1 (the minimum) 
XX=2.260 and 2.26 (the cell boundary) 
CNT=3 Three observations 

between 2.26 (see above) 
X¥x=2.420 and 2.42 (the next boundary) 
CNT=1 One observatio 

between 2.42 (eee aoe 
Xx=2.580 and 2.58 (the next boundary) 
CNT=3 Three observations 

between 2.58 (see ewe): 
Xx=2. 740 and 2.74 (the next boundary) 
CNT=1 One observation 

between 2.74 (see above) 
Xx=2.900 and 2.90 (the maximum) 


14. After the last cell boundary is displayed, the calcuia- 
tor will display "CMD" asking the user for the next command. 
Frequently, the user will wish to modify the histogram by 
changing the number of cells or the cell width. To reéecalcu- 
late the histogram cell frequencies without re-entering the 


data press: 


0 ALPHA AGAIN ALPHA 


If no further work with the histogram is desired, the user 


may view the order statistics simply by pressing run/stop. 
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mI LRODUCTION: 

Linear programming is an operations research technique 
normally associated with computerized data bases and the 
largest computers. Because of the complexity of the con- 
puter programs for linear programming and the large amount 
of data associated with most real world problems, calcula- 
tors have not been widely used for this application. With 
the increased memory capacity of the HPYICV, however, it is 
now possible to offer a calculator program which can solve 
interesting small scale linear programs. Of value primarily 
as an educational aid, this program will also be able to 
solve many small scales problems found at Army division, bri- 
gade and battalion level. The text by Hillier and Lieberman 
fRef. 8: pp. 16-66] is recommended as an introduction to 
the theory of linear programming as used by the program 
given in this appendix. Use of the program requires the 
user to formulate the linear programming problem; set upa 
Simplex tableau in standard form including adding slack, 


Surplus and artificial variables as required; and interpret 
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mie final tableau including the calculation of the values 
associated with the variables in the final basis. Using the 
Secleau form of the Simplex algorithm, the calculator per- 
forms both phase I (to obtain a feasible solution) and phase 
II (to obtain an optimal solution) to solve the linear pro- 
gramming problem. The calculator automatically determines 
the pivot column and pivot row for each pivot step. Infeas- 
ible and unbounded problems are automatically identified for 
the user by the program. There is no explicit handling of 


variables with upper bounds. 


PROGRAM DESCRIPTION: 

The program is written as a series of subroutines, each 
Seewhich performs a major step in the Simplex algorithm. To 
provide clarity to the user, alphabetic labels have been 
retained to identify the subroutines in lieu of faster and 
more memory efficient numeric labels. The alphabetic labels 
have not been retained for usSe aS program entry points and 
may be changed to numeric labels at the option of the user. 
meen program has two entry points, "LP" for running a new 
problems and "ALP" for reviewing data previously entered. 

Subroutine “"FINDQ" determines the pivot column by 
selecting the variable in the objective function with the 


most negative "price." If “PINDO" discovers at least one 
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negative value in the objective function, then the tableau 
column number associated with the most negative value will 
Bemstored in register 05. Upon return from “FINDOQ," the 
main routine tests register 05 to see if it contains a non- 
zero entry. If the entry is zero, it means that no further 
pivots will improve the value of the objective function, and 
the Simplex algorithm halts. If the program was in phase I 
(flag 11 clear) and the value of the phase I objective func- 
tion is reduced to zero, then a feasible solution has beén 
found and the program will automatically proceed +o phase ITI 
to discover an optimal solution. 

Subroutine "FINDP" determines which variable will leave 
the current basis by performing a minimum positive ratio 
test along the pivot column. In this way, the pivot row is 
determined. The row number of the pivot row is stored in 
register 06. Upon return from subroutine "FINDP," the main 
routine tests register 06 to see if it contains a non-zero 
entry. If the entry ls zero, it means that the problem is 
unkounded and the Simplex algorithm halts. Such an 
unbounded condition is most likely caused by an error in the 
problem formulation. 

Having determined the pivot column and the pivot fow, 


Subroutine "PIVOT" performs the actual Simplex pivot 
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operation. To speed calculation register 00 is used as a 
Pemmorary register to hold the reciprocal of the pivot 
element. Note that the pivot row is handled separately from 
the other rows in the tableau. lag 04 is used to provide 
the option of stopping calculation after every pivot. When 
this flag is set, the program will halt to allow the user to 
review the status of the tableau with the "ALP" function. 

Subroutine "CHECK" has two primary functions. First, it 
is used to verify that the designated basic variables are in 
row @limination form prior to the start of the Simplex 
algorithm. This means that the basic variable must have a 
coefficient of 1 in the row in which it is basic and zero's 
in all other rows. The second function of check is to pre- 
pare the objective function for phase I, if the initial 
basis contains artificial variables as indicated by one or 
more minus signs in the "JB" vector. 

Three subroutines are used to query the user for input 
data. Subroutine "READMN" queries the user for the number 
of constraints and decision variables in the problem and 
verifies the calculator memory is set to contain all the 
data necessary to solve the problem. Subroutine "READJB" 
queries the user for a column vector of pointers which indi- 


cate which variable is currently basic in each row. When 


de 
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Meecrang this vector Of pointers, the user indicates 
trtificial variables with 2 minus Sign. Subroutine "READA"™ 
queries the user for the values in the initial Simplex 
tableau including the slack and surplus variables and the 
right hand side and objective function. 

Several other service routines also are provided in this 
program. Memory size verification is done by subroutine 
Soman," which is called from within "READMN." Subroutine 
"In" is used to query the user for data entry and is called 
mewedll Of the data input routines. Subroutine "NXT" ini- 
tializes registers which contain frequently used quantities 
such as the the total size of the tableau for phase I and 
phase II. Subroutine "INIT" clears the calculator memory 
and sets flags and program constants appropriately for input 
of a new problem. Subroutine "SETL" establishes the loop 
counters used repeatedly within almost every other subrou- 
mene. Subroutine "“ERR1I" displays an appropriate error 


message if a data entry error is detected. 


SAMPLE PROBLEM: 
A division assistant G4 is planning an ammunition upload 


plan. There are four types of tank munitions to consider, 
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mec Luding: 


A = Discarding Sabot Rounds 

B = High Explosive Anti-Tank Rounds 
C = Phosphorous Munitions 

D = Machine Gun Ammunition 


Based on the Commander's guidance the assistant G4 is to 
consider the sabot rounds as twice aS important eae HEAT 
rounds, which in turn are themselves twice as important as a 
unit amount of phosphorous munitions and machine gun 


ammunition. His mission then, is to maximize: 
Z2= WA + 2B + C + D 


He ls, however, constrained by the following factors: 


“ 


ils There can be no more than 30 units of both sakot 
and HEAT munitions combined. 


2. There can be no more than 50 units of all types of 
ammunition combined. 


3. There must be at least 30 units of HEAT and 
phosphorous munitions combined. 


4. There must be at least 5 units of machine gun 
ammunition. 


These constraints may be expressed as: 


A+ 8B S30 
het wee. C+D Ss | 50 
Bae 230 

Daa 5 


Based on the Commander's guidance and the constraints listed 
above, formulate an optimum load plan. Fractional units are 


permitted. 


oie: 





SOLUTION: 

1. Before beginning with the calculator, the first step is 
to layout the tableau in standard form. This step and the 

last step of interpreting the final tableau requir2 working 
knowledge of linear programming as explained in Hillier and 
Liberman {[Ref. 8: pp. 16-66]. The standard form of the 


tableau is: 


1 2 3 4 cE 6 7 8 9 jOMmeiad 
A 3 C D H1 H2  S1  S2 Al A2- RHS 
1 1 30 
1 1 1 1 1 50 
1 1 -1 1 30 
4 -1 1 5 
Pa ea SSS sr sae ay 


In this tableau, H1 and H2 are surplus variables; $1 and S2 


are slack variables; and Al and A2 are artificial variables. 


2. The first step with the calculator is to set the size of 
the calculator's data memory. This program requires 20 reg- 
isters for temporary storage, 1 register for each tableau 
element, and 1 register for each row to hold the pointer to 
the baSic variable for that row. Thus, if M is the number 
of constraints and N is the number of variables including 


Slack, surplus and artificial variables, then the total data 
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storage requirement is: 


storage required = 21 + M + ((N + 1) x (M + 2)) 


AS mentioned in the program description, the "SIZE" subrou- 
tine will automatically verify that the user has allocated 
enough data storage to solve the problem. The length of the 
program is such that 177 data storage registers is the maxi- 
mum number of data storage registers that can be allocated. 
Thus, linear programs with 7 constraints and 15 decision 
variables can be solved with this program. For this 


example, press: 


= XEQ ALPHA Soe ALPHA 175 


3. Call for execution of the program with a new data set. 


Press: 


XEQ ALPHA Le ALPHA 


4. The calculator will respond with the prompt "NUM ROWS?" 
asking for the number of constraints in the linear progran 


formulation. In this example, press: 


G R/S 


5. The calculator will respond with the prompt "NUM COLS?" 


asking for the number of variables in the problem. The user 
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must count the number of slack, Surplus and artificial 


Variables in this total. In this example, press: 


LOS R7o 


6. The calculator will respond with the prompt "BASIC 1 ?" 
asking for the variable number of the variable which is 
basic in the first row. One of the major features of this 
program is that the basic variables need not be in the 
rightmost positions in the tableau. Thus, if a tableau were 
given in which some pivots had already been performed, the 
program could resume operation immediately. In ¢his 


example, press: 


Ea / Ss 


In a Similar fashion, the calculator will then query the 
user for the variable number of the variables which are 
basic in the remaining rows. 


For this example: 


sce Respond With 
Basic 2 ? 8 R75 
Basic 3°? 9 Cisne 7 
Basic 4 ? 10 CHS R/S 


Notice that because the basic variables in rows three and 
four are artificial variables, the variable number is 


entered as a negative number. This signals the calculator 
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that these variables must be driven from the basis in order 


#o reach an initial feasible solution. 


7. Next, the calculator will respond with "T1,1?" asking 
for the first element in the tableau. The usar should enter 
the numbers in the tableau using the digit entry keys and 
pressing run/stop after every entry. Notice that the right 
hand side and the objective function will be entered with 
the appropriate index in the tableau as shown in step 1 
above. The user must insure that the objective function is 
in standard form with the appropriate sign for each 


coefficient--in this example each coefficient is negative. 


8. After the last element in the tableau has been entered, 
the calculator will begin to automatically perform the Sin- 
plex algorithm. If the user wishes to stop the calculator 


after every pivot, he may at any time press; 


R/S SF O04 R/S 


If this flag is set, the calculator will stop and display 


the pivot number after every pivot is completed. 


9. When the Simplex algorithm can no longer improve the 
objective function, the calculator will stop and display the 


value of the objective function. In this example, the 
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calculator will stop after approximately three minutes and 


display: 


VALUE=110.000 


10. At this point, the user must use entry point "ALP" to 
determine the status of the final tableau. For this 


example, press: 


XEQ ALPHA ALP ALPHA 


Then by sequentially pressing the run/stop and clear arrow 
keys, the basic variables and tableau entries will be 


displayed. For example, in this problen: 


2¢ee Press see Meaning 

BASIC 1? CLX 2 Varaanule 2 is basic 
an che. £rest row. 

BeoLtc 2? CLX 1 Variable 1 is basic 


tmeethe Seconce LOW. 
etc. 


Then for the elements of the tableau: 


see Press see Meaning 
i, 1? CLX 0.000 Tableau entry 
aig 2 2 CLX 12000 Tableau entry 


11. After the calculator is finished, it remains for the 
user to interpret the final tableau. Again, the reference 
by Hillier and Lieberman [Ref. 8: pp. 16-66] is of primary 
value. In particular, the user must be able +o determine 


the value of the final decision variables based upon what 
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Weaeranles are in the basis, and what the final "right hand 
side" values are for each row. For this é¢xample, the final 


tableau is: 


» § @ bd #1 #2 $1 82 & a2 Rds 
-1  -1 1 -1 1 15 
1 1 1 1 -1 =-1 15 
1 1 1 -1 15 
1 -1 1 5 


eno. © om 2 2° |. 3S =-2 =-2  °&4°4¥°10 
Thus, the solution may be interpreted as 15 units each for 


Munitions A,B and C and 5 units for munition D. 
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me4iC SOURCE CODE: LINEAR PROGRAMMING 


THE FOLLOWING TABLE DESCRIBES THE KEY REGISTER AND FLAG 
ASSIGNMENTS MADE BY THIS PROGRAM: 


ROG — Monee nin ShoGioteRees HOLDS RECIPROCAL OF PIVOT 
ELEMENT IN SUBROUTINE PIVOT. 
ROQO1 = LOOP INDEX COUNTER 
ROQO2 = LOOP INDEX COUNTER 
RUS = TEMPORARY REGISTER. HOLDS MIN VALUE IN FINDOQ; 
MAX VALUE IN FINDP; AND IS AN INTERMEDIATE 
ADDRESS CALCULATION VALUE IN PIVOT. 
ROY = Aeon Be dane BOR SENDERECT ADDRESSING 
ROS = THE PIVOT COLUMN NUMBER 
RO6 = THE PIVOT ROW NUMBER 
RQ? = EFFECTIVE ZERO LEVEL 
RO8 = M = NUMBER OF ROWS 
RO9 = M PLUS 1 
R10 = M PLUS 
R11 = N = NUMBER OF VARIABLES 
Ri2 = N PLUS 1 
R13 = BASE REGISTER FOR THE LOCATION OF THE 
VECTOR JB WHICH CONTAINS POINTERS TO 
WHICH VARIABLE IS BASIC IN EACH ROW. 
R14 = ROW NUMBER OF THE OBJECTIVE FUNCTION; 
Serelert ells 1 OR MN PLUS 2 AS CETERULNED 
Diowo ED OR UPHASE 1. 
R15 = BASE REGISTER FOR THE LOCATION OF THE 
PHASE I OBJECTIVE FUNCTION, IF NEEDED. 
R16 = TEMPORARY REGISTER. 
Ri7 = NUMBER OF PIVOTS PERFORMED. 
Rilo = RESERVED FOR PULURE USE. 
R1i9-R177 = DATA STORAGE REGISTERS FOR ELEMENTS OF 
THE TABLEAU AND THE JB VECTOR. 
FLAGS 
FO1 - FO3 = SUBROUTINE EXECUTION FLAGS. BECAUSE THESE 
FLAGS ARE VISIBLE IN THE DISPLAY THEY CAN 
BE SET WHEN ENTERING A MAJOR SUBROUTINE AND 
CLEARED WHEN LEAVING ~~ GIVING THE USER A 
VISUAL INDICATION OF WHAT IS HAPPENING 
INSIDE THE CALCULATOR. 
FO1--SUBROUTINE FINDQ 
FO2--SUBROUTINE FINDP 
FOS =-SUBROUTINE PIVOT 
FO4 = WHEN SET, STOPS CALCULATOR AFTER EACH PIVOT. 
FO7 = USED AS TEMPORARY FLAG IN PIVOT AND CHECK ROUTINES. 
F10 = USED AS A TEMPORARY FLAG IN READ ROUTINES. 
Fil = WHEN SET, INDICATES PHASE II IS IN PROGRESS. 
B29 = SCONTROLS FORMAT OF DISPLAY SEPARATOR. 


483 END 
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PN TRODUCTIONS: 

The calculator subroutines described in this appendix 
perform functions which are frequently required by applica- 
tion programs and are therefore ideal candidates for use in 
a read only memory (ROM.) These routines were written espe- 
Cially to illustrate the differences between read only men- 
ory routines and routines designed for individual use via 
bax code or magnetic cards. These differences include more 
peecentiOn to entry and exit point options, an attempt to 
keep the size of the routines as compact as possible, and an 

f[oempt not to disturb the register stack if at all 
possible. 

These common subroutines are provided separately fron 
application programs because when more than one application 
program uses the routines, as is recommended for these pro- 
grams, the use of a separate block of common functions saves 
Space in the ROM overall. Also, by providing a convenient 
set of "macro" instructions, application prcgrams can be 


constructed more quickly and easily. Because these 
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Subroutines are used frequently, they have been individualiy 
optimized and tested to save memory space and execution 
time. By uSing «hese "macros" within an application pro- 
gram, the application programmer can be reasonably certain 
of their efficiency and reliability. 

Almost all user/calculator interface is handled by these 
routines. There is one set of Subroutines which assumes the 
user has a printer, and one set which does not. Printer 
instructions are preceeded in the listings shown in the 
appendix by an (PRT: label. When not using these routines 
on read only memory, the user will load one set or the other 
fore not both), as appropriate for his/her application. In 
so doing, the user with the printer gets full benefit from 
it while the user without the printer pays no penalty in 
execution time or memory space for the calculator's print 
mmetcructions. Also, to change from use of the printer to 
use of the calculator without it, the user needs only to 
read in the new common block--the application programs are 
retained in memory unchanged. The subroutines appear «he 
Same to any application program--giving the added benefit 
that any application program which uses them for input or 


output operations will automatically make good use of 
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the printer even if written by a programmer who did net 
axplicitly consider a printer requirement. 

When using these common subroutines, a discipline is 
enforced upcn the application program concerning use of the 
calculator memory registers. This saves the programmer 
from having +o plan his "register map" from scratch for each 
new program. Also, it insures compatability across differ- 
ent application programs for similar data objects such as 
Matrices and loop index counters. One of the most annoying 
protlems with read only memory programs available from the 
calculator manufacturer is this lack of cross program con- 
patability. Conflict in the use of memory registers is the 
rule, rather than the exception; and it is frequently impos- 
Sible to efficiently use more than one read only memory pro- 
gram as a subroutine within a user written program. A third 
reason why register asSignment standards are advantageous is 
that they make it easier for the user of the calculator to 
remember the key register assignments and, if necessary, 
Becall their contents during the execution cf an application 
progran. 

Another function performed by this set cf common subrou- 
tines 1s *0 simplify the use of indirect addre2ssing--a 


@eetical goal on the HPWicy. 
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Because the common subroutines listed in this appendix 
acre always called by application programs and never from the 
keyboard by the user, the typical user instructions are 
inappropriate. Instead, for the benefic of application 
programmers wishing to use the routines, the basic functions 
and structure of each are explained in subsequent sections 


of this appendix. 


SP a Be ee SS a oe 


Subroutine "IN® is used as a general input and output 
interface between the user and the calculator. This subrou- 
tine has three alternative entry points which when called 


affect functions as follows: 


ee mode (displays a_guestion mark fee 
IO--Output node ieee ate labeled data value) 
IX¥--Direct mode (input of value in x register) 


In particular, one entry point, "IN", may be called whenever 
an application program must query the user for a numeric 
input value. As such, it is a direct replacement for the 
PROMPT instruction crganic to the calculator, but automati- 
cally prompts, verifies and stores the received value using 
an indirect address contained in register 00. The printer 
version of the subroutine will automatically label and print 


the final, verified data value recorded. 
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Subroutine "IN" uses register 00 as a data location 
pointer and automatically increments this register so thar 
subseguent calls to the subroutine will result in sequential 
data manipulation. The application programmer must insure 
that register 00 contains a number equal to the storage reg- 
ister number prior to calling the subroutine. For example, 
if register 00 contains 17, "IN" will store the data in reg- 
ister 17. One of the major advantages of this routine is 
that the same subroutine may be used to verify and/or change 
the data previously recorded. Thus, separate edit routines 
are usually unnecessary. Pressing the R/S key without 
touching any other key leaves the value stored unchanged. 
Pressing "1" and *#"* and then "R/S" adds one «*o the current 
stored value, and so on for other function keys. Entering a 
new string of digits results in that new value being stored. 

An additional feature of this subroutine is the "verify" 
mode indicated by flag 05. Flag 05 is reserved for this 
purpose and is set "on" by a call to subroutine "VRI'"-- 
another of the subroutines in this common set. The verify 
mode is intended for use by a novice or i user who 
wishes to verify every data value as it is keyed into the 
calculator. The advantage is increased accuracy and 


confidence in the result. 





Subroutine "D2" 
Subroutine "D2" is used to set up the ind2x register for 
a program loop. This subroutine has two alternative entry 


points which when called, increment different index regis- 


ters as follows: 


D2--Establishes register 02 as the index 
Di--Establishes register 01 as the index 
This subroutine is intended for use with the "ISG" loop 


structure which has the effect most like that of a FORTRAN 


"DO LOOP." For example, to execute a loop 20 times: 


rtetoe ¢ 
ObIi@=s e« 
rt 
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00 


00 
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eo © FIN 
©? OM 


The advantage of this form of loop structure is that 
register 00 may be used within the loop for address calcuia- 
tions. The first time the loop is addressed the integer 
portion of the number in register 00 will be 1, the second 
memge it will be 2 and so on. There is no need to truncats 
the fractional portion of the number because the HP4YIC 
ignores the fractional component of a number when 


calculating a register address. Use of index registers for 


address calculation makes indirect addressing practical. 
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Registers 01 and 02 should be reserved for use as index 
registers by the application programmer. In most cases 


these two registers should prove sufficient. 


Subroutine "VR" is used as a general purpose calculator 
initialization routine. This subroutine has «three alterna- 
tive entry points which vary the amount of initialization 
performed as follows: 


VR--Sets the verify mode on, and the following; 
WR--Suppresses the audio tones, and the following: 
WS--Clears all. memory ; 

Sets the display for integers, 

Assigns statistical regisers 

Sets "zero" level for équalit tes 


5 qd 
5 : CLUNg ¢ Nn 
Sets base addr2ss for indirect address 


a 
DniNyee 

Mire ehe printer version of this initialization routine, 
the subroutine prints a banner (usually the program name) 
which has been stored in the alpha register prior to calling 
the initialization subroutine. 

Peet lag 13 1S set pEr#or to "calling the initialization 
routine, then the calling pregram must have placed the nuna- 
ber of data registers required to execute the program in the 
mMebegister prior to calling the initialization routine. [In 
this case, a check will automatically be performed using 
subroutine "SZ" described below. 

It is recommended that all application programs provide 


an alternate entry point which bypasses the initialization 
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step. Then if a data error is encountered, the user nay 
review the data entered into the calculator by simply press- 
ing the return key once after every prompt. This procedure 
works because subroutine "IN" recalls the stored value prior 
momprOompting the user. When the user presses the clear key, 
the alphabetic prompt is removed and the existing data value 


revealed. 


Subroutine "sz" 

Subroutine "SZ" is used to test if sufficient numbers of 
data registers are available +o run an application progran. 
This subroutine may be either called directly or as part of 


mie initialization routine described above. 


euoLoutine "ER" 

Subroutine "ER" is called whenever the application pro- 
gram encounters an error--usually in the input data. A 
prompt is displayed and an audio tone sounded, provided flag 
26 (the audio enable flag) has not been cleared by the 


initialization routine described in paragraph D. 


Subroutine "SORT" 
This subroutine is included to illustrate the use of a 
stack register table in the program comments, but it also 


represents a useful utility routine. The sorting algorithm 





Pacedees tre spell sore [Ret. 6: pp. 84-95] which gives 
reasonably fast sorting times with a very small progran 
size. All conventions such as base register in ROY and nun- 
ber of data points in R15 are assumed by this subroutine. A 
complete list of all such register assignments is listed at 


the end of the program listing. 


These two small routines provide a useful capability to 
store and recall up to three integers between 0 and 999 in 
one data register. This means that if you ar2 manipulating 
a spread sheet of small, positive integers you can store the 
Same data in one third the space. Of course, run «ime is 
degraded (about 20 seconds for every 100 data references.) 
To store a value, assuming the base register has been 


defined, just press: 


value ENTER| point-number XEQ "PUT 


To recall a value, simply press: 


point-number CE OmeuGr 
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DESIGN METHODOLOGY: 

This appendix discusses the design methcdology used dur- 
ing construction of both the cross compiler program and the 
command processor, which is an IBM EXEC II program which 
provides an interactive programming environment for users of 
the systen. 

Blazie*s compiler for the HP65 calculator [Ref. 9] rep- 
resents one of the first attempts to provide a compiler for 
Seeeulator programs. Both Carvalho {Ref. 10: pp. 25-29] 
and McNeal [{Ref. 11: pp. 148-178] have published BASIC lan- 
guage programs which cross compile HP4Y1CV instructions on a 
microcomputer for output to a line printer which can rrint 
acceptable bar code. While these referenced programs pro- 
vided valuable insights into the problem, especially into 
the special characteristics of the HP41CV instruction set, 
none was 2xactly suited to the needs of this study. Pecause 
the Versatec plotter at the Naval Postgraduate School could 
be easily used only by FORTRAN programs, FORTRAN seemed the 


computer language of choice for this project. Both programs 
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were written with limited objectives and neither would have 
easily suppcrted the extensions desired. Extensions planned 


pereeanplementation included: 

e An extended instruction set. 

e In code comments. 

e Extensive error checking. 

e Compatability with the Emulator. 
Meo ynthetic Instructions [ Ref. 1]. 


Sernscruction macro's, 

Having decided to code an original cross compiler, a 
design methodology which would capitalize on the advantages 
of FORTRAN was planned. FORTRAN'sS major deficiency for use 
mconstriucting a compiler of any type is its lack of 
alpha-numeric string handling capabilities. Rather than 
Struggle with the lack of string functions, it was decided 
to code the necessary string functions as separate subrou- 
tines. This decision was reinforced countless times 
BEaroughout the process of writing the compiler. The string 
function subroutines have been used not only in the cross 
compiler, but in many other FORTRAN programs since they were 


Originally written. In fact, many persons who have no 


ee 





interest in the HP41CV cross compiler may find the Set of 
string functions listed in this thesis to be a valuable set 
of utility routines to be used to augment FORTRAN. The gen- 
eral convention used througout the string function subrou- 
tines is that an alphabetic string may be represented as a 
vector of two byte integer variables used to store the char- 
acters and a single four byte integer variabl= ased to 
store the length of the string. 

One of the major advantages of the cross compiler is its 
ability to handle comments integrated within «he HP4YI1CV 
source code. This feature is critical to the clarification 
of the legical structure of the HP4Y1ICV programs. Because 
the parenthesis is not a valid HP41CV character, it was cho- 
sen as the comment indicator character. A comment may occur 
beginning at the first column on an input line or anywhere 
after an HP4Y1ICV instruction. The comment must follow the 
instruction because everything after the comment mark out to 
the end of the input line 1s considered part of the comment. 

The control the user has over the output listing is also 
one of the advantages of the cross compiler. When twe con- 
ment indicator marks are placed in positions one and two of 
the input line, the compiler will force a page eject when 


Printing the output listing. In addition, the user can vary 
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the number of output lines per page and cause useful tanners 
to be placed adjacent to program labels for ease of 
recognition. 

Altogether there are twenty-four subroutines and a main 
program which constitute the cross compiler. The source 
code for each of these routines is provided in the second 
section of this appendix. Each subroutine begins with a 
Statement concerning its function and construction. <aAccord- 
ingly, no general description of éach subroutine will be 
repeated here. However, subroutine COMP deserves special 
attention, for it is the master lexicographic analyzer for 
the compiler and would also interface the user with the enu- 
lator. Its function is to receive a single line of HP4ICV 
source code and identify it. COMP considers all HP4I1CV 
instructions to be of one of three types. The firs* cat- 
egory are the single byte instructions with no operands that 
can be compiled by a simple table look up. COMP has been 
constructed so that the instruction set can be extended at 
any time simply by increasing the size of this table. In 
this way abbreviated or altered command names could be eas- 


ily used. The second category of instructions are th 


( 


multi-byte instructions which require a table lookup and the 


translation of one or more operands, including possibly an 


1 





Magirect 1nStruction indicator. The table examined by the 
compiler is the same as for the category one instructions, 
and a cod2 is given in the table which indicates to the 
compiler the number of operands which are required with each 
instruction. A syntax check 1s then made in subroutines 
IONE and ITWO to insure that the number and characteris- 
tics of each operand are appropriate for the given instruc- 
tion. One of the major advantages of the use of the cross 
compiler is the syntax and error checking that is performed 
during the compilation process. The third type of instruc- 
tion represents the exceptional instructions that are so 
difficult +o compile that they require separate subroutines 
for efficient compilation. These instructions include stor- 
age and recall of data, program labels and program flow con- 
trol statements such as goto and execute. 

In order to provide an efficient programming command 
system for the compiler that would minimize the need to know 
technical details about the operation of the compiler, an 
IBM EXEC II program was written. This program no* only 
interfaces the user to the compiler, but it also provides on 
line user instructions as to how to use the system. 

Included in this command processor is a command menu which 


gives the format and short description for each command. 
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Another command, help, provides more detailed information 
about each command. When a novice user first enters the 

exec, or types the name of the exec followed by a question 
mark, then he receives a four page narative description of 
what the system is, how it works, and what actions he must 


take to write a successful HPY1CV application program. 
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